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We update our numerical investigation of topological structures around static quarks in pure gauge QCD by 
results of the first runs including dynamical quarks. Simulations were performed on an 8 x 4 lattice, with SU(3) 
Wilson action, with 3 flavours of quarks of equal mass, both in the confinement and deconfinement phase. In the 
confinement phase we observe indications for the existence of a flux tube between a static quark and antiquark, 
flux-tube breaking for large separations, and local correlation between the topological charge density and chiral 
condensate. In the deconfinement phase almost all configurations turn out to be topologically trivial. 



1. The LATTICE 93 contribution of our group 
[1] presented results of the study of topological 
properties of QCD vacuum configurations with 
static colour sources in pure SU(3) lattice gauge 
theory. The present work extends that investiga- 
tion to SU(3) with dynamical quarks. The main 
results from pure gauge QCD simulations can be 
summarized as follows [1,2]: 

- Configurations with nonvanishing topological 
charges appear very seldom in the deconfinement 
phase. 

- The topological charge density is lowered in the 
vicinity of external colour sources in both phases 
of pure gauge QCD. 

- The suppression of the topological charge den- 
sity appears in the whole flux tube between the 
static quark and antiquark in the confinement 
phase, while in the deconfinement phase the effect 
is concentrated around the (anti)quarks only. 

It is of direct interest to verify if the above pic- 
ture remains unchanged after taking dynamical 
quark degrees of freedom into account. In gen- 
eral one expects that most qualitative features of 
the confinement mechanism are the same in QCD 
without and with dynamical quarks. However, 
in a theory with fermions the topological suscep- 
tibility is proportional to a power of the quark 
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mass [3]. Therefore, one can reasonably expect 
that the observed topological effects will depend 
on the presence of dynamical quarks and their 
masses. 

2. Our calculations were performed on an 
8 3 x 4 lattice with (anti)periodic boundary con- 
ditions using the Metropolis algorithm. We sim- 
ulated QCD with SU(3) gauge group, with stan- 
dard Wilson action for gauge fields, and with 3 
flavours of Kogut-Susskind quarks of equal mass, 
ma = 0.1. We performed runs for two values 
of j3, one for the confinement (/? = 5.2) and the 
other for the deconfinement phase (/? = 5.4). In 
the former simulation observables of interest were 
measured on 2000 configurations separated by 50 
sweeps, in the latter run only 300 configurations 
have been measured until now. The parameters 
of our runs were chosen in such a way so that the 
corresponding physical length scales were approx- 
imately the same as in our pure gauge simulations 
[1,2], thus enabling direct comparison of observed 
physical effects. 

Topological properties of lattice configurations 
were probed by two local lattice operators of 
topological charge density: 
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Top. charge probabilities, plaquette definition (beta=5.2) 
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Figure 1. Probabilities for various values of the 
topological charge at /? = 5.2. 



Figure 2. Squared topological charge density 
around a static quark at /? = 5.2. 
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(the "plaquette" (P) and the "hypercube" (17) 
operator [4]). 

These operators are known to differ from the 
continuum topological charge density by mul- 
tiplicative renormalization factors substantially 
different from unity at the values of /? used [5]. 
These factors were removed using the cooling 
method, in particular a slightly modified variant 
of Hoek et al. [6]. 

The local distribution of the topological charge 
density around single quarks and around a me- 
son with quark-antiquark separation d was ob- 
tained by calculating the correlation functions 
(L(0)q 2 (r)) and (L(O)L^ (d)q 2 (r)), respectively, 
and normalizing them to the corresponding vac- 
uum expectation values. L(r) as usual denotes 
the Polyakov loop describing the propagation of 
a static colour charge. 

3. The results of the simulation of full QCD at 
/? = 5.2 are very similar to those of pure-gauge 
calculations at /? = 5.6 [2]. After a few cooling 
steps topological charges of lattice configurations 



measured using the plaquette and hypercube op- 
erators coincide and reach nearly integer values. 
The probability distribution of the charges Q is 
shown in Fig. 1. It is approximately Gaussian as 
in the pure-gauge case (cf. Fig. 2 of Ref . [2]). 

The topological charge density (squared) is 
suppressed in the vicinity of a quark (Fig. 2). The 
suppression again occurs in the whole flux tube 
between the quark and antiquark (Fig. 3a, b). 
However, Fig. 3c reveals an essential difference 
between pure-gauge configurations and those con- 
taining dynamical quarks. At the largest QQ sep- 
aration, d = 4, the suppression is weaker in the 
middle of the tube than in the pure-gauge sim- 
ulation. This we consider as a sign of flux-tube 
breaking due to the creation of a virtual dynam- 
ical qq pair. A similar effect was reported earlier 
in Ref. [7]. 

Our data also show correlations between the 
topological charge and the chiral condensate. The 
confinement configurations with higher topologi- 
cal charges tend to have higher values of the chiral 
condensate. The correlation appears also locally 
which is clearly seen from comparing Fig. 2 with 
Fig. 4, depicting the local chiral condensate dis- 
tribution around a static quark. 

The deconfinement simulations (/? = 5.4) are 
still in progress. Most configurations belong 
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(a) 5 cooling steps, d = 2 



Figure 4. Local distribution of the chiral conden- 
sate around a static quark at /? = 5.2. 
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(b) 5 cooling steps, d = 3 



to the topologically trivial sector, and non-zero 
charges appear even more rarely than in the pure- 
gauge case. 

4. Our results from first runs with dynamical 
quarks indicate similar behaviour of topological 
quantities in pure-gauge QCD and full QCD, and 
hence a similarity in the corresponding confine- 
ment mechanisms. The quark-mass value used, 
however, was rather high and our conclusions 
have to be supported by simulations at lower mass 
values. 




(c) 5 cooling steps, d = 4 



Figure 3. Squared topological charge density 
around a static quark-antiquark pair at /? = 5.2. 
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